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We present a new analysis of double charmonium production in e+e^ annihilation. The observation 
of the processes e^e~ J/iprjc, J/ipXcO, and J/ip rjc{2S) is confirmed using a dataset more than 
three times larger than that of Belle's previous report, and no evidence for the process e^e~ — > 
J/ip J/xp is found. We perform an angular analysis for J/ipr^c production and set an upper limit on 
the production of J/ip J/'tp. Processes of the type e^e~ — > tp{2S){cc)rcs have been observed for the 
first time; their rates are found to be comparable to those of e^e~ J/ip{cc)rcs processes. 

PACS numbers: 13.66.Bc,12.38.Bx,14.40.Gx 

The large rate for processes of the type e^e^ — )■ J/ijjrjc and J/ip (cc)non-rcs reported by Belle Q remains 
unexplained. Following the publication of this result, the cross-section for e+e^ J/iprjc via e^e^ 
annihilation into a single virtual photon was calculated using non-relativistic QCD (NRQCD) to be 
~ 2fb 2], which is at least an order of magnitude smaller than the measured value. Several hyjDotheses 
have been suggested in order to resolve this discrepancy. In particular, the authors of Ref. 3] have 
proposed that processes proceeding via two virtual photons may be important. Other authors (Jl suggest 
that since the dominant mechanism for charmonium production in e+e^ annihilation is expected to be 
the color-singlet process e^e~ — s- ccgg, the final states observed by Belle contain a charmonium state 
and a Af ~ 3 GeV/c^ glueball. Such glueball states are predicted by lattice QCD and can have masses 
around 3GeV/c^. Possible glueball contributions to the Xco signal are also discussed in Ref. Q. 

The previous Belle analysis was performed with a data sample of 45 fb^^. The process e+e~ J/ipVc 
was inferred from the rjc peak in the mass spectrum of the system recoiling against the reconstructed J/ip 
in inclusive e+e^ — > J/ip X events. In this paper we report an extended analysis of the e+e^ — > J/ip (cc)ios 
process to check the above hypotheses and provide extra information that might be useful to resolve the 
puzzle. This study is performed using a data sample of 140 fb~^ collected at the T(45) resonance and 
15 fb"^ at an energy 60MeV below the T(4S'). The data were collected with the Belle detector |^ at the 



KEKB asymmetric energy e+e" storage rings 

The analysis procedure is described in detail in Ref. ^ . For J/i/i reconstruction we combine oppositely 
charged tracks that are both positively identified either as muons or electrons. For J/?/; e+e~, the 
invariant mass calculation includes the four-momentum of photons detected within 50 mrad of the 
directions, as a partial correction for final state radiation and bremsstrahlung energy loss. The J /ip — *■ 
signal region is defined by a mass window \Mi+i- — Mji^\ < 30MeV/c^ (w 2.5(7 a/). QED processes 
are significantly suppressed by the requirement that the total charged multiplicity (A^ch) in the event be 
A^ch > 4. The contribution from J/tp mesons in BB events is removed by requiring the center-of-mass 
(CM) momentum P*j^^, to be greater than 2.0 GeV/c. A mass-constrained fit is then performed to improve 
the P*j/^ resolution and the recoil mass M^-c 



(Ecu - ^P*jh '^^ calculated, where 



J/ip 



the J /'ip CM energy after the mass constraint. ip{2S) is reconstructed via its decay to J/i/; tt+tt and the 
■0(25') signal window is defined as \Mj/^^+^-^ - M{'ip{2S))\ < lOMeV/c^ (w StTAf)- 

The Mrccoii(>//V') spectrum for the data is presented in Fig. Q clear peaks around the nominal ijc and 




recoil 



(J/¥) 



(GeV/c^) 



FIG. 1: The mass of the system recoiling against the reconstructed J/ip in inclusive e'^e J/ip X events. The 
curves are described in the text. 

Xco masses are evident; another significant peak around ^ 3.63GeV/c^ is identified as the ric{2S). The 
authors of Ref. 3] estimated that the two-photon-mediated process e~^e~ — > J/ipJ/ip has a significant 
cross-section and suggested that the observed e+e^ J/ipr]c signal in Q might also include double J/ip 
events, thereby producing an infiated cross-section measurement. Since e~^e~ annihilation to J/^ J/tp 
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via a single virtual photon is forbidden by charge conjugation symmetry, it was ignored in our previous 
analysis. To allow for a possible contribution from the exchange of two virtual photons, we fit the spectrum 
in Fig. ^ including all of the known narrow charmonium states. In this fit, the mass positions for the 
Vc, XcO and ric{2S) are free parameters; those for the J/ip, Xcij Xc2 and i-!{2S) are fixed at their nominal 
values. The expected line-shapes for these peaks are determined from a Monte Carlo (MC) simulation 
as described in our previous paper Q, the background is parameterized by a second order polynomial 
function, and only the region below the open charm threshold (A/iccoii < 3.7GeV/c^) is included in the 
fit. The fit results are listed in Tabled The yields for 77c, XcOi and rjc{2S) have statistical significances 
between 3.8 and 10.7. The significance of each signal is defined as ln(£o / C,max) , where Cq and 

'Cmax denote the likelihoods with the corresponding signal yield fixed at zero and at the best-fit value, 
respectively. The fit returns negative yields for the J ftp and i/'(2S'); the Xci and Xc2 yields are found to 
be consistent with zero. A fit with all these contributions fixed at zero is shown as a solid line in Fig. ^ 
the difference in the rjc, XcO and rjd'^S) yields compared to the default fit is small, and is included in the 
systematic errors. The dashed line in the figure corresponds to the case where the contributions of the 
J/i': Xci, Xc2 and ■0(25) are set at their 90% confidence level upper limit values. The dotted line is the 
background function. 

TABLE I: Summary of the signal yields (N), charmonium masses (M), significances, and cross-sections {a Bom x 
S>2((cc)ros)) for e^e~ — > J/V' (cc)rcs; S>2 denotes the branching fraction for final states with more than two 
charged tracks. 

(cc),es N M[GeV/c2] Signif. crBo,„xS>2 [fb] 

235 ± 26 2.972 ±0.007 10.7 25.6 ± 2.8 ± 3.4 

J/V- -14 ±20 fixed — <9.1at90%CL 

XcO 89 ±24 3.407 ±0.011 3.8 6.4 ±1.7 ±1.0 

Xci±Xc2 10 ±27 fixed — <5.3at90%CL 

77^(25) 164 ±30 3.630 ± 0.008 6.0 16.5 ± 3.0 ± 2.4 

V'(2S) -26 ± 29 fixed — < 13.3 at 90% CL 



Given the arguments in Ref. '3'|, it is important to check for any momentum scale bias that may 
shift the recoil mass values and confuse the interpretation of peaks in the A/iccoii spectrum. We use 
e+e^ ijj{2S)^, 0(25) — > J/^tt+tt" events to calibrate and verify the recoil mass scale. Events with 
a reconstructed ^{28) and with no other charged tracks form a pure e+e^ ip{2S)^ sample with less 
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than 1% background as estimated using the ipiiS) sideband region. We use the ip{2S) momentum to 
calculate the square of the mass of the recoiling system; the resulting spectrum is shown in Fig. [3 The 




FIG. 2: The square of the mass of the system recoihng against the reconstructed '0(25') — > J/ip tt'^tt' , in events 
with charged multiplicity equal to 4. Points with error bars show the data, the solid line shows the result of 
the fit described in the text. The hatched histogram (scarcely visible) shows the spectrum in the scaled ip{2S) 
sidebands. 

scaled ip{2S) sideband is also shown. We perform a fit to the M'^^^^-^^['tl:{2S)) spectrum, using Monte 
Carlo simulation to determine the expected signal shape; second order QED corrections, which produce 
a higher M^^^^^^ tail, are taken into account. The peak position is left free in the fit and the non-?/'(2S') 
background is ignored. The fit finds the shift in the data with respect to the MC function to be consistent 
with zero (AAfj^^^^jj — 0.010 ± 0.009 GeV^/c^). From this result we conclude that the J/tp recoil mass is 
shifted by not more than 3MeV/c^ in the region M^ccoi\{J /i^) ^ 3GeV/c^. 

As an additional cross-check we fully reconstruct double charmonium events. The rjc is reconstructed as 
— > TT+TT ) or 2{K'^K ) combinations within a window of ±50 MeV/c^ around the nominal 
7]c mass. In events with Nc\i = 6 we find 3 events with J /ij) rjc combinations in a ±100 MeV window around 
the CM energy 3 cr). No events are seen in the rjc sideband region (100 < AI (K s Ktt / 2{K K)) — M,,^ < 
350 MeV); a fit to the mass distribution gives an rjc signal significance of 4.1(t. Based on the rjc yield in 
the Afrccoii(^/'0) distribution, we expect 2.6 ±0.8 fully reconstructed events, consistent with the observed 
signal. Thus we conclude that the peak in M^ccoii{ J / ip) is dominated by rjc production. We also search 
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for fully reconstructed double J /ij) candidates in events with iVch = 4. No J/ if} J candidates are found 
in a window of ±100 MeV around the CM energy. 

Based on the calibration of the Af recoil ('^/V') scale, the result of the fit to the M-^^coiiiJ 1 4') distribution 
and the full reconstruction cross-check, we confirm our published observation of the process e'^e^ — > 
J/il^rjc and rule out the suggestion of Ref. that a significant fraction of the inferred J/'!/'?7c signal 
might be due to J/-0 J/i/; events. 

The reconstruction efficiencies for the J/iprjc, J/4'XcO, and J/ipr]{2S) final states strongly depend on 
^prod, the production angle of the J/^ in the CM frame with respect to the beam axis, and the helicity 
angle 6'hGi, defined as the angle between the decay £^ direction and the boost direction of the CM frame 
in the J/tp rest frame. We therefore perform an angular analysis for these modes before computing cross- 
sections. We fit the M,-ccoii{J / 'fp) distributions in bins of |cos(0pi.od)| and |cos(0hci)|, and correct the yield 
for the reconstruction efficiencies determined bin-by-bin from the MC. The results are plotted in Fig.|2| 
together with fits to functions A{1 + acos^O) (solid lines). We also perform simultaneous fits to the 
production and helicity angle distributions for each of the (cc)ics states, assuming J/tp (cc)i.cs production 
via a single virtual photon and angular momentum conservation, thus setting apiod = cthci- The values 
of the parameter a from the separate fits to | cos(0hci)| and | cos(0pi-od)|, and from the simultaneous fits, 
are listed in Table Hll 



TABLE II: The a parameters obtained from fits to the production and helicity angle distributions for e^e 







Separate fits 




Simultaneous fits 


(CC) res 


Qlprod 




Qhcl 


Ohcl = aprod 


XcO 


-1.7 ±0.5 




n 5+° '^ 
-0 7+° '' 

"•'-0.5 



U.O_Q 7 


0.93«j? 

-i.oilgj^ 

0.87t°:«^ 



The angular distributions for the J/tprjc and J / ijj rjc{'2,S) peaks are consistent with the expectations for 
production of these final states via a single virtual photon, Ofprod — ahoi — +1 1^- There is no evidence 
for the sharp rise in cross-section for large | cos(0prod)| expected for J /ijj J /ip production via two virtual 
photons 01 . The prediction for a spin-0 glueball contribution (e+e^ — > J/ipGo) to the J/^-q^ peak, 
ttprod — ckhoi — —0.87 Q, is also disfavored. 
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FIG. 3: Distributions of cosines of the production (left) and J/tp helicity angles (right) for e+e" J/tj) rjc (top 
row), e^e" — > J /ip XcO (middle row) e^e~ — > J rj{2S) (bottom row). The solid lines are results of the individual 
fits; the dotted lines are the simultaneous fit results. 

The process e^e^ —+--)'*—!• J/ijjXcO can proceed via both S- and D-wave amplitudes, and predictions 
for the resulting angular distributions are therefore model dependent. Our results disfavor the NRQCD 
expectation aprod = cthci — 0.25 and arc more consistent with S -wave production, where ckprod — 

ahoi = -1- 

To calculate the cross-sections for the processes e+e^ — > J/iprjc, J/tpXcO, we fix the pro- 

duction and helicity angle distributions in the MC to 1 -1- cos^ 9 for J / ip rjc{rjc{2S)) ^ and to 1 — cos^ 9 
for JjipXca. The statistical errors in the a parameters for the angular distributions are translated into 
uncertainties in the efficiency determination and included in the systematic error. To set a conservative 
upper limit for e^e~ ~f J/ipJ/ip, J/'<PXci{2)i '^/V' '0(2»S'), we use assumptions for the production and 
helicity angle distributions that correspond to the lowest detection efficiency. Note that for J/ -0 J/ip and 
J/ipip{2S) the assumed angular distributions lead to lower efficiencies than those that follow from the 
predictions of Ref . jSj . 
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To reduce the model dependence of our results due to the effect of initial state radiation (ISR) , whose 
form-factor dependence on of the virtual photon is unknown, we calculate cross-sections in the Born 
approximation. We first calculate the fraction of events in the signal Mi-cconiJ / ip) distributions that are 
accompanied by an ISR photon of an energy smaller than a cutoff i?cutoff using MC. We then correct the 
cross-sections calculated for that cutoff value using a factor that yields the Born cross-section 8]. The 
final result (ctbo™ ~ 0.70 • <yfuii) is then independent of the choice of the cutoff energy provided -Ecutoff 
satisfies Me ^ £^cutoff <C Ecm- As in Ref. because of selection criteria we present our result in terms 
of the product of the cross-section and the branching fraction of the recoil charmonium state into more 
than 2 charged tracks: a x S>2, where S>2((cc)rcs) = 'S((cc)ros — ^ > 2 charged). The cross-sections are 
given m Table D 

We perform a similar study with reconstructed ip{2S) J/i/; tt+tt" decays to search for e+e" 
ip {'2 S){cc) res processes. The recoil mass spectrum for the data is presented in Fig.0] peaks corresponding 
to the T]c, XcOi and •qd'iS) can be seen. The fit to the Afi.ccoii('0(25')) distribution is identical to the 




2.2 2.6 3 3.4 3.8 



Mrecoii(¥(2S)) (GeV/c') 

FIG. 4: The mass of the system recoiling against the reconstructed tpi'iS) in inclusive e'^e~ ip{2S)X events. 
The curves are described in the text. 

-Mrocoii(>//V') fit; but due to the limited sample in this case, the masses of the established charmonium 
states are fixed to their nominal values; the rjc{2S) mass is fixed to 3.630 GeV/c^ as found from the 
Afrocoii(>//^) fit. The signal yields are listed in Table Hill Significances for the individual 77^, XcOj f^nd 
ric{2S) peaks are in the range 3 ^ 4cr; the significance for e+e~ — > i^i^S) (cc)ios, where (cc)i.os is a sum 
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over rjc, XcO: E^nd r]c{2S), is estimated to be 5.3 cr. The significance is calculated as ln(£o / -Cmax) , 

where Cq and £max denote the likehhoods of the fit with all signal yields fixed at zero and at the best-fit 
value, respectively. In Fig. 01 the result of a fit with only rjc, XcO sj^id ric(2S) contributions included is 
shown as a solid line; the dashed line shows the case where the J/'0, Xci: Xc2, and ip{2S) contributions 
are set at their 90% confidence level upper limit values. The dotted line is the background function. 

To estimate the efficiency we assume the ipi^S) production and helicity angle distributions to be the 
same as those for the corresponding J/ip (cc)ros final states. Finally, the calculated products of the Born 
cross-section and the branching fraction of the recoiling charmonium state into two or more charged 
tracks (ct x S>o, where K>o((cc)rcs) = B{{cc)i-cs charged)) are presented in Table ITTll 

TABLE III: Summary of the signal yields (A'^), significances, and cross-sections (aso™ x Byoiiccjrcs)) for e'^e" — + 
V'(2S) (cc)rcs; S>o denotes the branching fraction for final states containing charged tracks. 



(cc)res N Signif. CTBorn X S>0 [fb] 

r^c 36.7 ±10.4 4.2 16.3 ± 4.6 ± 3.9 

J/il^ 6.9 ±8.9 — < 16.9 at 90% CL 

Xco 35.4 ±10.7 3.5 12.5 ±3.8 ±3.1 

Xci + Xc2 6.6 ±8.0 — < 8.6 at 90% CL 

T]c{2S) 36.0 ±11.4 3.4 16.0 ±5.1 ±3.8 

'tp{2S) -8.3 ±8.5 — < 5.2 at 90% CL 



The systematic error is dominated by the fitting systematics of the signal yields: 10% for J/tp (cc)ios 
and 14% for 'ip{2S){cc)rcs- To estimate this contribution we vary the parameterizations of the signal 
(intrinsic widths of charmonium states, form- factor dependence on Q^) and background in the fit to the 
-^recoil spectra. Another large contribution is due to the reconstruction efficiency dependence on the 
angular distributions (7% for J/ip (cc)rcs and 15% for ■ip{2S){cc),.cs)- In the MC, the angular parameters 
a are varied within the statistical errors of our angular analysis for e+e^ — s- J/^p (cc)rcs and in the 
full range (—1 < aprod(Q!hoi) < 1) for e+e^ '>P{2S) (cc)i-cs to estimate the uncertainty in efficiencies. 
Other contributions for e+e~ — > J/ip{cc)rcs (V'(2S')(cc)i.os) come from the multiplicity cut (3%(2%)), track 
reconstruction efficiency (3% (5%)) and lepton identification (3% (3%)). 

In summary, using a larger data set we confirm our published observation of e~^e~ J/iprjc^ J/ipXcO 
and J / T]c{2S) and find no evidence for the process e+e" J / ip J / if) . We have calculated the cross- 
sections for e^e^ J/^prjc, J/i^XcOj and J/^pric{2S) with better statistical accuracy and reduced sys- 
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tematic errors and set an upper limit for (T(e+e~ — > J j t\} J j ■\\}) x B{J/-4' — > > 2 charged) of 9.1 fb at 
the 90% CL. Although this limit is not inconsistent with the prediction for the J/t/j J/i/; rate given in 
Ref. 01 , the suggestion that a large fraction of the inferred J / ip rj^ signal consists of J/?/' J /ip events is 
ruled out. We have measured the production and helicity angle distributions for e+e^ ^ J /i}) rjc, J/4' XcO: 
and J/ipr]c{2S); the distributions are consistent with expectations for these states, and disfavor a spin-0 
glueball contribution to the rjc peak. We observe '0(25')(cc)ros production for the first time, and find 
that the production rates for these final states are of the same magnitude as the corresponding rates for 

J/'0 (cc)rcs- 
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